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ABSTRACT 

While the metaUicity excess observed in the central stars of planetary systems 
is confirmed by all recent observations, the reason for this excess is still a subject 
of debate : is it primordial, is it the result of accretion or both? The basic 
argument against an accretion origin is related to the mass of the outer convective 
zones which varies by more than one order of magnitude among the considered 
stars while the observed overabundances of metals are similar. We show here 
that in previous discussions a fundamental process was forgotten : thermohaline 
convection induced by the inverse /^-gradient. "Metallic fingers" may be created 
which dilute the accreted matter inside the star. Introducing this effect may 
reconcile the overabundances expected in case of accretion with the observations 
in stars of different masses. 

Subject headings: exoplanets; stars : abundances; stars: accretion; diffusion; 
double-diffusive convection ; thermohaline convection 

1. Introduction 

The metallicity distribution of the central stars of planetary systems compared with 
other stars of the same spectral types clearly show a metallicity excess of a factor two on the 
average, while the individual [Fe/H] values lie between -0.3 and -1-0.4 (Santos et al 2003 and 
references therein). Although this excess is generally discussed in terms of iron abundances 
only (for most observers [Fe/H] is a synonymous of "metallicity"), spectroscopic observations 
show that the abundances of other heavy elements beginning with carbon are also enhanced 
(Santos et al 2001) 
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The proposed possible explanations for this behavior are of two kinds: 

1) accretion hypothesis : this assumes that the star had normal abundances at its 
formation, but that it accreted metal-rich matter during the phases of planetary formation 
(Gonzalez 1998) 

2) primordial hypothesis : this assumes that the stellar gas out of which the planetary 
system formed was already metal-enhanced. A corollary in this case would be that such a 
metal enrichment is necessary for planets to form around stars 

Other explanations have been proposed which mix the two previous scenarios. 

The strongest argument invoked against the accretion hypothesis is related to the mass 
of the outer convective zone which varies by more than one order of magnitude for stars 
between 1.0 and 1.4 M0, while the metallicity excess is the same. On the other hand the 
primordial hypothesis has to deal with the fact that some of the host stars have a sub-solar 
metallicity, which is an argument against the assumption that metallicity excess is needed 
for planetary formation. 

Diffusion of accreted metals has been discussed by several authors (Pinsonneault et al 
2001, Murray and Chaboyer 2002). They show that gravitational settling is not rapid enough 
to produce a decrease of the metal abundances and lead to the observed values. They also 
discuss rotational and other mixing processes in relation with the observed lithium abun- 
dances. As lithium is detroyed by thermonuclear reactions at a relatively small temperature 
(3.10^K) inside stars, it is often used as a test for mixing processes. However such a com- 
parison between the metal diffusion and the lithium depletion is highly model dependent as 
it assumes a specific variation of the mixing effect with depth. 

Here we show that the argument against accretion does not hold as soon as we take into 
account the convective instability induced by inverse /x-gradients, or thermohaline convec- 
tion (also called "double-diffusive convection"). If freshly accreted metals accumulate in a 
small convective zone on top of radiative layers, they will not stay there but will be diluted 
downwards in "metallic fingers" similar to the "salt fingers" observed in the ocean, so that 
their overabundance will rapidly decrease with time. 

We discuss how this process should take place, with the conclusion that the final metallic 
abundances should not depend on the depth of the standard convective zone as usually 
assumed. However the exact amount of metals which may remain in the stellar outer layers 
depends on parameters like the size and depth of the metallic fingers, which cannot be 
precisely constrained in the framework of our present knowledge. We only show that a final 
overabundance by a factor two as observed can be obtained with plausible values of the 
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unknown parameters. Two ways can be explored to go further in this study : numerical 
simulations of metallic fingers and asteroseismology of overmetallic stars. 



Thermohaline convection is a wellknown process in oceanography : warm salted layers 
on the top of cool unsalted ones rapidly diffuse downwards even in the presence of stabilizing 
temperature gradients. When a blob is displaced downwards, it wants to go further down 
due to its overweight compared to the surroundings, but at the same time the fact that 
it is hotter contradicts this tendancy. When the salt gradient is large compared to the 
thermal gradient, salted water normally mixes down until the two effects compensate. Then 
thermohahne convection begins. While the medium is marginally stable, salted blobs fall 
down like fingers while unsalted matter goes up around. This process is commonly known 
as "salt fingers" (Stern 1960, Kato 1966, Veronis 1965, Turner 1973, Turner and Veronis 
2000, Gargett and Ruddick 2003). The reason why the medium is still unstable is due to the 
different diffusivities of heat and salt. A warm salted blob falling down in cool fresh water 
sees its temperature decrease before the salt has time to diffuse out : the blob goes on falling 
due to its weight until it mixes with the surroundings. 

The salt finger instability can occur with any two components which have different 
diffusivities if there is an unstable gradient of the slower diffusive component and a stable 
gradient of the faster diffusive component. For this reason it is also referred to as "double 
diffusive convection". In oceanography, the fastest diffusive component is conventionally 
referred to a T while the slowest one is S even when the components are different. Indeed 
the effect is generally studied in the laboratory using water with a mixture of salt and sugar. 
In this case sugar is the slowest component as it diffuses about three times slower than salt. 
The results can well be visualized, pictured and studied (e.g. Wells 2001). The fingers have 
also been studied by 2D and 3D numerical simulation (e.g. Piacsek and Toomre 1980, Shen 
and Veronis 1997, Yoshida and Nagashima 2003) 

The condition for the salt fingers to develop is related to the density variations in- 
duced by temperature and salinity perturbations. Two important characteristic numbers 
are defined : 

• the density anomaly ratio 



2. 



Thermohaline convection 



2.1. 



The salted water case 



Rp = aVT/pVS 



(1) 
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where a — — (^^)5 p and ^ — (^||)t,p while VT" and V-S" are the average temperature and 
sahnity gradients in the considered zone 

• the so-called "Lewis number" 

T = Ksji^T = TtIts (2) 

where ^^s and ^^T are the saline and thermal diffusivities while rs and tt are the saline and 
thermal diffusion time scales. 

The density gradient is unstable and overturns into dynamical convection for Rp < 1 
while the salt fingers grow for Rp > 1. On the other hand they cannot form if Rp is larger 
than the ratio of the thermal to saline diffusivities t"^ as in this case the salinity difference 
between the blobs and the surroundings is not large enough to overcome buoyancy (Huppert 
and Manins 1973, Gough and Toomre 1982, Kunze 2003). 

Salt fingers can grow if the following condition is satisfied : 

l<Rp< r'^ (3) 

In the ocean, r is typically 0.01 while it is 1/3 for a salt-sugar mixture. We will see below 
that in solar-type stars where T is the temperature while S is the mean molecular weight 
this ratio is about 10~^° if k,s is the molecular (or "microscopic") diffusion coefficient but 
it can go up by many orders of magnitude when the shear fiow instabilities which induce 
mixing between the edges of the fingers and the surroundings are taken into account. 



2.2. The stellar case 

Thermohaline convection may occur in stellar radiative zones when a layer with a larger 
mean molecular weight sits on top of layers with smaller ones (Kato 1966, Spiegel 1969, 
Ulrich 1972, Kippenhahn et al 1980). In this case = dln///dlnP plays the role of the 
salinity gradient while the difference Vad ~ V (where Va,^ and V are the usual adiabatic and 
local (radiative) gradients dlnT/dlnP) plays the role of the temperature gradient. When Vad 
is smaller than Vrad, the temperature gradient is unstable against convection (Schwarszchild 
criterion) which corresponds to warm water below cool water in oceanography. In the op- 
posite case the temperature gradient is stable but the medium can become convectively 
unstable if (Ledoux criterion): 



Verit = JV^ + Vad - V < 



(4) 
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where (f) — {d hi p/d In fj) and 5 — {d In p/d In T) When this situation occurs, convection 
first takes place on a dynamical time scale and the enriched matter mixes down with 
the surroundings until Vcrit vanishes. Then marginal stability is achieved and thermohaline 
convection may begin as a "secular process" , namely on a thermal time scale (short compared 
to the stellar lifetime!). 

Such an effect has previously been studied for stars with a helium-rich accreted layer 
(Kippenhahn et al 1980). It was also invoked for helium-rich stars in which helium is sup- 
posed to accumulate due to diffusion in a stellar wind (as proposed by Vauclair 1975) and 
for roAp stars in case some helium accumulation occurs (Vauclair et al 1991). Similar com- 
putations have been done in the case of accretion of matter on white dwarfs, in relation with 
novae explosions (e.g. Marks and Sarna 1998). 

As we will show below, if hydrogen poor matter is accreted on the top of a main- 
sequence type star with normal abundances, it creates an inverse /^-gradient which may lead 
to thermohaline convection. Comparing the stellar case with the water case, we can guess 
that metallic fingers will form if the following condition is verified : 

l<|^'^°--^'|<x-' (5) 

with r = Dfi/DT = tt/t/j, where Dt and Z^^ are the thermal and molecular diffusion coeffi- 
cients while Tt and are the corresponding time scales. This condition is similar to condition 
(3) for the stellar case. In the following, we will neglect the deviations from perfect gas law, 
so that 4> = 6 = 1. 

In the next section we show computations of the time scales and orders of magnitude 
of this process. 



3. The fate of accreted metals in solar type stars 

We study the case of main-sequence solar type stars which would have accreted hydrogen- 
poor material at the beginning of their lifetime. We assume, for simplicity, that the accretion 
occurred in a very short time scale compared to stellar evolution. We then study the fate of 
the accreted metals and chose the examples of 1.1 M© and 1.3 Mq stars. 

As the chemical composition of the accreted matter is not known, we assume that all 
elements are accreted with solar relative abundances except hydrogen and helium which are 
assumed completely absent. We will see that our general conclusion would not be changed 
if the relative abundances were modified : the time scales would only be slightly different. 
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After accretion, the metal dilution occurs in two phases. Phase one : rapid convection 
takes place in a dynamical time scale until it reaches the marginal equilibrium obtained when 
Vcrit vanishes (Equation 4 ). Phase two : thermohaline mixing begins and dilutes the metal 
excess until condition (5) is no longer satisfied. 

In the following, we first discuss the //-gradients induced by the original metal excess 
in the convective zone. We compute the situation at the end of phase one, according to the 
value of the accreted mass. We discuss the case of the same accretion occurring for the two 
different stellar masses. Finally we study the time scales of thermohaline mixing and the 
remaining overabundances in the convective zone. 



The main parameters of the two stellar models we have used as an example are given 
in table 1. They correspond to 1.1 Mq and 1.3 Mq, both at an age of about 1.5 Gyr. We 
can see that the convective zone is 7 times more massive for the 1.1 M0 model than for the 
1.3 M0 model, which does not change sensitively during main sequence evolution. We thus 
expect that metal accretion leads to an original overabundance 7 times larger in the 1.3 Mq 
model than in the 1.1 Mq model. 

For completely ionised hydrogen and helium, the mean molecular weight may be ob- 
tained with the following expression : 



Here A represents an average mass of metals, (^) a relative abundance with respect to 
hydrogen, and x an averaged number of particules (ion and electrons) associated with the 
metals. In the following we take (^) = 0.1 and (^)o = 1.4.10"^ for the original abundances 
of hehum and metals (Grevesse and Sauval 1998). We treat the metals as an average element 
to which we attribute the mass A = 20. We assume hydrogen and helium completely ionized 
and we neglect the term x{^). In these conditions we obtain /xq = 0.6 for the value of the 
mean molecular weight before accretion. If metals are added, /j, is modified by : 



3.1. 



Dynamical convection and marginal equilibrium 




(6) 



2 + 3(f)+a;(f) 



Aji 




(7) 



Let us now write : 




M 



M 
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and define ctj as the initial value of a, obtained just after the accretion process (assumed 
rapid compared to the other time scales). 



The variations of (^) and fi are then related by : 



M 

A/. = 9(a-l)(-)o (9) 



We have computed, as a function of the initial overabundance ratio ai, the depth at 
which metal-enriched material has been diluted when it reaches the marginal equilibrium 
phase, and the actual overabundance ratio a in the convective zone. This has been obtained 
by integrating in all cases the mass of metals diluted by this process. 

Suppose, for example, that the metal excess ratio in a 1.1 M© star is a = 2 after dilution 
by dynamical convection. This corresponds to an original overabundance factor of 2.05. The 
same accretion mass would lead to an original metal excess ratio ai = 14.3 in a 1.3 Mq star. 
Our computations show that such an original metal excess is diluted in a "transition zone" 
down to ttz = 7.5.10^'^cm, leading to an overabundance after dilution a = 8.5. Figure 1 
displays the abundance profiles obtained in this case in the 1.3 Mq star. The parameters at 
the bottom of the transition zones as we have defined above are shown in table 1 for the two 
models. Figure 2 and 3 display, for the two models, the values of and the corresponding 
q; as a function of the width of the transition zone. 

Thermohaline mixing begins after this dilution process. 
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Table 1. Model data' 



model 


r 






T 


Pr 


Pr 




n/i(yr) 


1.1 (CZ) 


5.2el0 


0.75 


3.7e31 


2.0e6 


.128 


3.5el3 






1.1 (TZ) 


4.9el0 


0.71 


4.6e31 


2.3e6 


.195 


6.2el3 


15.8 


4200 


1.3 (CZ) 


S.OelO 


0.84 


5.1e30 


9.0e5 


0.007 


1.2el2 






1.3 (TZ) 


7.5el0 


0.79 


8.2e30 


1.2e6 


.013 


2.2el2 


12.6 


800 



'^the values are given at the bottom of the convective zone (CZ) and at the 
bottom of the "transition zone" (TZ) defined as the region of metal dilution 
after dynamical convection, for the example discussed in the text (metal 
excess a = 2 in the 1.1 Mq model) ; the metal abundances decrease with 
depth inside this region until it reaches its normal value at ttz (cf. Figure 

1)- 
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Fig. 1. — Profiles of tlie metallic overabundance factors in the 1.3 Mq model before (ctj) 
and after {a) dilution by dynamical convection for the example discussed in the text, as a 
function of the fractional radius inside the star. Thermohaline convection begins after this 
phase. 
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3.2. Thermohaline mixing, metallic fingers and time scales 

The study of thermohaline mixing in stars is far from trivial. Detailed comparisons 
of numerical simulations and laboratory experiments in the water case have recently been 
published (Gargett and Ruddick 2003) but the stellar case may be different as mixing then 
occurs in a compressible stratified fluid. 

In the following we use the formalism proposed by Kippenhahn et al 1980 (hereafter 
KRT) who discussed the time scale of thermohaline mixing in stars in a simplified way and 
applied it to the case of helium rich layers standing on the top of hydrogen rich ones. 

In first approximation they treated thermohaline mixing as a diffusion theory. The 
whole picture may be described as follows: blobs of metal enriched matter begin to fall down 
from the convective zone and exchange heat and heavy elements with their surroundings; 
chemicals diffuse more slowly than heat, so that the blobs go on falling down until they are 
completely disrupted, thereby creating finger shapes. The most efficient process for element 
diffusion out of the blobs is the shear flow instabifities at the edge of the flngers : as the 
falling down matter undergoes friction with matter going up, turbulence occurs and mixes 
part of the fingers with their surroundings, on a horizontal length scale which is a fraction 
e of the horizontal size L of the blobs. The blobs disappear when they have travelled down 
on a distance W long enough for this mixing to disrupt them completely. 

KRT define a diffusion coefficient Dth as the product of the blobs vertical velocity Vfj, by 
their size L. They evaluate as : 

V, = ,^ \. i^) (10) 

where Dji represents the ji difference between the inside and the outside of the blobs. The 
diffusion coefficient D^h is then obtained as : 

Hp dlnn 
|V„d - V| T*f^H dr 

where Hp is the pressure scale height and t^j^ the local thermal time scale given by : 

16acr3 ^^^^ 

where the parameters have their usual meanings. We can see that vanishes in this diffusion 
coefficient which does not depend on the size of the blobs. 

In this simple picture, Dth is not exactly the local mixing coefficient as it involves 
the whole blobs while local mixing involves only the edges of the blobs. The local mixing 
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Fig. 2. — This figure shows, for the 1.1 model, the initial metallic excess ctj = (^)i/ (^)o 
and the metallic excess a obtained after dilution by dynamical convection, as a function of the 
depth {tcz — r) at which they are diluted. The a curve is obtained with Wcrit — (Equation 
4). For example, an original overabundance of 30 will be diluted inside a transition region 
of thickness 1.05.10^°cm, below the standard convective zone. An overabundance of 15 will 
remain after this process. Then thcrmohaline convection will go on reducing this value. The 
curves close to the origin are displayed in a zoom. 
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coefRcient becomes correctly represented by Dth only at the bottom of the fingers, where 
the blobs disappear. For this reason, KRT define the thermohaline diffusion time scale as : 
^th = / Dth where W is the depth of the fingers (it also represents the size of the transition 
zone in which ix varies by A/x, between its two extreme values). This diffusion time scale 
may be written : 

^'""-^^-I^^ ^ ^^^^ 

Note that if the mixing scale at the finger edges increases linearly while the blobs fall down, 
the coefficient e is nothing else than the ratio of the horizontal and vertical length scales of 
the fingers, namely L/W. 

We give in table 1 the thermohaline diffusion time scales corresponding to the example 
given in the previous section, with a metallic excess of a factor 2 in 1.1 M©. Computations 
done with a metal excess of 1.5 instead of 2 give similar results for the thermohaline diffusion 
time scales. They are always very short compared to the stellar lifetime. If nothing could 
stop these metallic fingers from growing, the convective zone would be emptied of its metallic 
matter short after it is accreted. 

We must now take into account the second condition for the formation of the fingers 
(Equation 5), namely : 

1^1 <- ^ (") 

At the end of the dynamical convection phase, the | ^ | ratio is one. Then, during 
thermohaline convection, the local /x-gradient decreases so that the | ^ | ratio increases. 
Thermohaline convection should stop when this ratio reaches the value of the Dt/D^ ratio. 
The remaining overabundances in the stellar outer layers strongly depend on the efficiency 
of the turbulent mixing at the edge of the fingers, which is non linear and not really known. 
We can only give here a tentative discussion of this process. 

The thermal diffusion coefficient Dt may be evaluated as the square of a relevant length 
scale (e.g. the size of the falling blobs) divided by the local thermal time scale t^^ : 

D - ^ - ^^"^^^ (15) 

Below the convective zone in the 1.1 Mq , Dt is typically of 3.0.10®cm^.s~^ while it is 
8.10^cm^.s~^ in the 1.3 Mq model. Meanwhile the microscopic element diffusion coefficient 
is of order one cm^.s~^ : if no turbulence occured, the metallic fingers would extend far down 
and deplete the overmetallic material until an extremely small //-gradient would be reached 
: in this case no metallic excess would be left in the convective zone. 
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If we take into account the shear flow instabihties which result from the motion of the 
blobs, has to be replaced by the local turbulent diffusion coefficient which is not exactly 
Dth but t^Dth where e represents the L/W ratio. Using equation (11) for Dt^, we find from 
condition (14) that the fingers disappear when W is of the same order as L, which makes 
sense but is not helpful to evaluate the amount of heavy elements which remains in the outer 
layers at the end of the whole process : the KRT approximations are too rough to be useful 
in this respect. 

Another approach consists in estimating the value of which would be needed to 
account for the observations. Our aim is that a mctallicity excess of order two remains 
in the convective zone after these processes occurred. From Equation (9) we find that it 
corresponds to A// ~ 0.0126 or, with = .6 as computed in these stellar layers, A/x//x ~ 0.02. 
If we suppose that the mixing region extends over one pressure scale height (which is by not 
means proved !) we find that the Dt/D^ ratio should be of order 10, which seems reasonable 
compared to the water case. Such a value would be obtained with a turbulent diffusion 
coefficient at the edge of the fingers of order lO'^ to 10^ cm^.s"^. 

More computations are needed for a better understanding of the shear flow turbulence 
induced by the flngers. Furthermore the growing of the metallic flngers may depend on the 
other processes at work in the star : rotation-induced turbulence, internal waves, etc. (cf. 
similar problems in salted water : Gargett and Ruddick 2003). 

In any case the computations presented here show that the metallic matter accreted 
onto a star does not stay inside the standard convective zone : it flrst turns over due to 
dynamical convection and then goes on diffiising due thermohalinc convection. The observed 
overabundances in planetary systems host stars can be obtained with realistic values of the 
unknown parameters. 

4. Conclusion and future prospects 

The important conclusion of this paper is that the strongest argument against accretion 
as an explanation of the metallic excess observed in host stars planetary systems host stars 
has to be revised : if hydrogen poor matter is accreted in the early phases of stellar evolution, 
during planet formation, the metal excess in the convective zone leads to a destabilizing /x- 
gradient which induces metal dilution, first by dynamical convection, then by thermohalinc 
mixing, on a time scale much shorter than the stellar lifetime (typically 1000 yrs). The 
remaining overabundances in the convective zones depend on the physical conditions inside 
the star. We do not expect, in case of accretion, to end up with metallic overabundances 
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increasing with decreasing convective mass as obtained in standard models. Furthermore, 
the region mixed by the thermohahne process is locahsed below the convective zones and 
should not go down to the lithium destruction layers. 

In the present study, we have supposed that the accreted material had a relative so- 
lar abundance, except for hydrogen and helium which were assumed completely absent. 
Changing these relative abundances would modify the relation between the ji value and the 
metallicity (Equation 9). The general conclusions would still hold : the numerical values 
and time scales would only be changed accordingly. 

We did not discuss the amount of matter needed to explain the observed abundances 
by accretion and how this accretion could have proceed. In the example of section 3.1, a 
metal excess of 2.05 in 1.1 Mq star would be obtained with 145 of all metals included, 
corresponding to about 11 of iron, which is a large value. A metal excess of 1.5 would 
be obtained with 5 of iron or 66 of all metals. The difficulty of explaining such a 
large accretion from protoplanetary discs represents a second argument against the accretion 
hypothesis, which we do not discuss here : this is a different subject, out of the scope of the 
present paper. 

Detailed spectroscopic observations and precise abundance determinations will be of 
great interest to test the accretion processes against the primordial scenario for exoplanets 
hosts stars. Asteroseismic studies of these stars will also help deriving whether they are 
overmetallic down to the center or only in their outer layers and give an important hint for 
a better understanding of the formation of planetary systems. 

I thank the referee N. Murray who gave very useful comments and critics on a first 
version of this paper and helped improving the discussion. 
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